J/ip Production in ^/sNN = 200 GeV Cu+Cu Collisions 
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Yields for J/i/> production in Cu+Cu collisions at ^/snJi = 200 GoV have boon measured over the 
rapidity range \y\ < 2.2 and compared with results in p-\-p and Au+Au collisions at the same energy. 
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The Cu+Cu data offer greatly improved precision over existing Au+Au data for J/-)/) production in 
collisions with small to intermediate numbers of participants, in the range where the Quark Gluon 
Plasma transition threshold is predicted to lie. Cold nuclear matter estimates based on ad hoc fits to 
d+Au data describe the Cu+Cu data up to A^'part ~ 50, corresponding to a Bjorken energy density 
of at least 1.5 GeV/fm^ 

PACS numbers: 25.75.Dw,12.38.Mh,21.65.Qr,25.75.Nq 



High-energy heavy-ion colhsions provide the opportu- 
nity to study strongly interacting matter at very high en- 
ergy densities where Quantum Chromodynamics (QCD) 
predicts a transition from normal nuclear matter to a de- 
confined system of quarks and gluons, the Quark-Gluon 
Plasma (QGP) Jj. At the Relativistic Heavy Ion Collider 
(RHIC) the energy density in central Au+Au collisions 
is well in excess of the critical energy density expected 
for this transition 

Over the past twenty years, there has been intense the- 
oretical and experimental work on J /ip production. First 
predicted by Matsui and Satz Q , suppression of quarko- 
nia production in ultra-relativistic heavy ion collisions 
was expected to be an unambiguous signature for the 
formation of a QGP. It is now recognized that in order 
to interpret J/ij] production as a QGP probe one has to 
consider cold nuclear matter effects such as initial state 
energy loss and shadowing , as well as charm quark 
energy loss [6 1 , co- mover interactions Q , corrections for 
feed-down from higher mass charmonium states, and sec- 
ondary production mechanisms, such as recombination of 
initially uncorrelated cc pairs [8| . 

Experiment NA50 reported suppression of J /if) pro- 
duction in Pb+Pb collisions at ^Jsnn = 17.3 GeV [9j 
that exceeds expectations based on their measurements 
of cold nuclear matter effects m. p + A collisions [lo| . 
NA60 observed similar behavior in In+In collisions at 
the same energy [ll|. The PHENIX experiment [12] at 
RHIC has characterized effects of the nuclear medium on 
J/^ production at ^/SNN = 200 GeV. The basic invari- 
ant yield reference is obtained bom p+p data 13, 3, 
Cold nuclear matter effects are studied using rf+Au 
data [3, Cold and hot nuclear matter effects are 

studied for large numbers of participants (A^part) using 
Au+Au data [l3,[l3|, and for smaller A^part using Cu+Cu 
data, the subject of this paper. The results are presented 
as a nuclear modification factor, i?AA, the ratio of the 
yield in heavy ion collisions to the yield in p+p collisions 
scaled by the number of binary nucleon-nucleon collisions 
(-^cou), which is appropriate for point-like processes. 

Lattice QCD calculations [l^ indicate that the thresh- 
old energy density for QGP formation is of order 1 
GeV/fm'^. At ^'snn = 200 GeV this is expected to occur 
below A^part = 100 in a region where Au+Au data 
have limited statistical and systematic precision [l8| . 
High statistics measurements with the intermediate sized 
system Cu+Cu provide crucial information in that im- 
portant region. 



In this Letter we present results obtained by PHENIX 
during the 2005 RHIC run on the production of J/ip 
in Cu+Cu collisions at y/sJTW = 200 GeV. J/ip invari- 
ant yields were studied via J/^ e+e~ decays mea- 
sured at midrapidity with the central arm spectrometers 
(|y| < 0.35, A(l) = 2 X 90°), and J/^p fi+p' decays 
measured at forward rapidity with the two muon arm 
spectrometers (1.2 < \y\ < 2.2, Ap = 360°). Event 
centrality and the location of the collision vertex along 
the beam axis (z^tx) are measured with two Beam-Beam 
Counters (BBC) located at 3.0 < \ri\ <3.9. A Glauber 
model and a simulation of the BBC response was used 
to determine A^pait and A'coii and their systematic uncer- 
tainties for different collision centrality ranges [21] . 

Data were recorded using lepton triggers in coincidence 
with a minimum bias trigger which required a coincidence 
between the BBC detectors and a valid z^tx- After apply- 
ing a cut of \zytx\ < 30cm and quality assurance criteria, 
the data correspond to a sampled luminosity of about 2.1 
nb^^ (1.3 nb~^) in the e+e" {fX^ analysis. 

Electron detection at midrapidity used the Drift 
Chambers for momentum measurement, the Pad Cham- 
bers for pattern recognition and track location, and the 
Ring Imaging Cherenkov (RICH) detector plus Electro- 
magnetic Calorimeter (EMCal) for electron identifica- 
tion. Charged particle tracks were matched with a RICH 
ring and an EMCal hit to select electron candidates by 
requiring at least two RICH phototube hits inside an an- 
nulus around the projected ring center, ring quality cuts, 
track/cluster position matching cuts at the EMCal, and 
a cut on the ratio of EMCal energy to track momentum, 
E/p - {E/p) > -2cr. 

The J/ij} — > e^e^ trigger required one signal above a 
certain energy threshold in the EMCal and a matching 
RICH hit. Two energy thresholds were used during the 
run, 1.1 GeV and 0.8 GeV, yielding average J/ip trigger 
efficiencies of ~ 65% and 82%, respectively. The J/ip ^ 
e^e~ signal extraction method was very similar to the 
method used in the recent Au+Au [31 and p + p [lB| 
analyses. The like sign invariant mass spectrum was sub- 
tracted from the unlike sign spectrum. The remaining 
yield in the J/i/' mass region (2.9 < Minv < 3.3GeV/c^) 
was corrected for pairs lost to the radiative tail and pairs 
added by the continuum signal under the peak [3| • The 
total J/tfj count in the e+e" channel was ~ 2,050. The 
signal to background ratio (S/B) was « 1(6) for the most 
central(peripheral) collisions. 

Muon detection at forward and backward rapidities 
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used the muon arms, consisting of cathode strip tracking 
chambers in a magnetic field (MuTr) and larocci tube 
planes interleaved with thick steel absorbers (MuID). 
Muon candidates were identified by penetration to the 
last MuID gap, and their momenta were measured by 
their bend through the MuTr. 

The dimuon trigger required two candidate tracks to 
penetrate the MuID, point back to the event vertex, and 
pass an opening angle cut {6 > 19°). The dimuon com- 
binatorial background was estimated using the product of 
the like sign counts, 2V N++ ■ N^^ , and was subtracted 
from the unlike sign spectra. The residual background 
(notably from the open charm pairs and Drell-Yan) was 
evaluated using an exponential form. The J/^ —^ fJ-'^ fJ-^ 
signal was estimated by direct counting of the remain- 
ing pairs above the exponential fit in the mass range 
2.6 < Minv < 3.6GeV/c^ and also by using two fits with 
different parameterizations (single and double Gaussian) 
of the J/V' line shapes, as described in [l^, The 
average of the results gave the signal count and the vari- 
ation gave the systematic error. The total J/tp yield was 
« 9,000. The S/B was w 0.3(1.0) for the most cen- 
tral(peripheral) collisions. 

The J/f/; invariant yield in the appropriate centrality, 
rapidity and transverse momentum bin is given by : 



J/V 



1 



N 



27rpT dp^dy 27rpT ^cvt AyApx^e ' 



(1) 



with Bu the branching ratio for J/V' Nj/^ the 

number of observed J /ip] N^^t the number of events; Ay 
the rapidity range; Apx the transverse momentum range, 
and Ae the acceptance and efficiency correction (includ- 
ing trigger efficiency). 

The determination of Ae is done with a full GEANT 
simulation. The method is described in more detail 
in [l^. Ae decreases with the collision centrality due 
to overlapping hits in the RICH and the EMCal in the 
central arm, and in the MuTr for the forward arms, lead- 
ing to an increasing fraction of misreconstructed tracks 
in higher multiplicity events. This effect is evaluated by 
embedding simulated single J/ip events in real events. 
The efficiency loss in the most central collisions is 3% 
for dielectron measurements and 20% (16%) for dimuon 
measurements at positive (negative) rapidity. 

Systematic uncertainties in the measured J/ijj invari- 
ant yield depend on J /ij) rapidity and transverse momen- 
tum as well as on event centrality. Systematic uncer- 
tainties are grouped into three categories: point to point 
uncorrelated uncertainties (type A), which can move the 
points independently of each other, point to point corre- 
lated uncertainties (type B), which can move the points 
coherently, though not necessarily by the same amount, 
and global systematic uncertainties (type C). In all plots 
point to point uncorrelated systematic uncertainties and 
statistical uncertainties are quadratically summed and 
represented by vertical bars, point to point correlated 



systematic uncertainties are represented by boxes, and 
global systematic uncertainties (if any) are quoted. 

TABLE L Systematic error sources, values and types for -Raa 
vs Apart in the two rapidity intervals. Where a range is given, 
it is from peripheral to central collisions. 



source 


\y\ < 0.3 


y\ G [1.2,2.2] 


type 


signal extraction 


6 % 


5-6 % 


A 


detector + trigger efficiency 


1.4-5 % 


3 % 


B 


run by run variation 


5 % 


2 % 


B 


input y + pt distributions 


2 % 


3 % 


B 


A^coll 


14-11 % 


14-11 % 


B 
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FIG. 1: (color online) J/t/j yield vs px at mid (left) and for- 
ward (right) rapidity for different Cu-|-Cu centrality bins and 
for p + p Uncertainties are described in the text. 

Systematic uncertainties of type A and B for Raa vs 
A'part are summarized in Table U Some uncertainties 
in the invariant yield, such as that on the acceptance, 
cancel out for Raa and are not shown. Global systematic 
uncertainties for Raa vs Apart include the p+p J/V' yield 
uncertainty and some p + p systematic errors that do not 
cancel when forming Raa- 

Results for the two muon arms agree within uncertain- 
ties and are combined where appropriate. Fig. [1] shows 
the J/il) yield vs pt for different Cu-f Cu centrality classes 
at mid and forward rapidity. As was done previously for 
the Au-I-Au case [3], the mean square transverse mo- 
mentum, (px), was calculated numerically from the data 
for pt < 5 GeV/c. The Cu+Cu data are plotted vs 
Apart and compared with the corresponding values from 
Au-hAu [la], d+ku [la] and p+p collisions in Fig. [2 
Within uncertainties, the data for Cu-|-Cu and Au-|-Au 
agree where they overlap in A^part, and the {p\) for the 
Cu-|-Cu data seems independent of A^part- 

The Raa values vs px and rapidity are shown in Fig. [3] 
for the 0-20% most central Cu-f Cu collisions. We see 
similar behavior for mid and forward rapidity, and there 



5 



> 

0) 

(5 
m 

V 



^ 2 




|y| < 0.35 
Cu+Cu 
O P+P 
□ d+Au 
■ Au+Au 



|y|£ [1.2,2.2] 
• Cu+Cu 
O P+P 
□ d+Au 
■ Au+Au 



10 



10' 



N, 



10 



10^ 



part 



part 



FIG. 2: (color online) The (px) A^'part for J/ip production 
in Cu+Cu, P + P Ts^, d+Au lS\ and Au+Au [3 collisions 
at mid (left) and forward (right) rapidity. 



J 2 
t 

1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 
Q 



|y| < 0.35, syst =±16% 



Centrality [0,20]% 




J/>|/- 
J/i|/- 



■ e e 



syst 



Centrality [0,20]% 



2 3 4 5 -2 

[GeV/c] 



-10 12 

Rapidity, y 



FIG. 3: (color online) Raa vs pt (left) and y (right) for J/t/j 
production in the most central Cu+Cu collisions. 



appears to be no dependence in all centrality classes. 
The RMS width of the rapidity distribution (evaluated 
directly from the data) is identical, within ^ 2 — 3% un- 
certainties, in P + P collisions and in all centrality classes 
for Cu+Cu collisions. 

Figures Hl^a) and [Ifb) show similar behavior within 
uncertainties for Raa in Cu+Cu and Au+Au 
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lisions at comparable values of iVpart- Theoretical calcu- 
lations [1^ including only modified initial parton distri- 
bution functions and an added J — N breakup cross 
section were fitted in [3] to d+Au J/ij} Raa data. The 
EKS98 [13, m and nDSg ^ shadowing models were 
used. The fit was made simultaneously to all rapidities 
by optimizing the breakup cross section. While consis- 
tent with the low statistics d+Au data [l^, this method 
leads to a model dependence of the CNM effects, since the 
rapidity shape is determined entirely by the shadowing 
model. In an attempt to reduce this model dependence, 
we used a data-driven ad hoc model to parameterize the 
d+Au data [ll|. The ad hoc model uses EKS98 (method 
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FIG. 4; (color online) (a,b) Raa vs A^'part for J/ip production 
in Cu+Cu and Au+Au [la ] collisions. The curves are predic- 
tions from ad hoc fits to d+Au data [3| and are discussed 
in the text, (c) Ratios of the measured -Raa values to the 
predicted cold nuclear matter i?AA. The dashed lines show 
the 1 a uncertainties from the d+Au fits. 



1) and nDSg (method 2) shadowing parameterizations 
for the relative rapidity dependence within the fitted ra- 
pidity ranges, but the breakup cross section is replaced 
with a quantity, which we call /, that is optimized sep- 
arately for y=0 and |y|=1.7. The fits using method 1 
yielded fdAu = 2.3±f;g mb at y=0 and 3.9±J;2 mb at 
|y|=1.7. The method 2 fits yielded fdAu = 0.9±}:| mb at 
y~Q and 3.3±};2 mb at |?/| = 1.7. The resulting separate 
parameterizations of the d+Au data vs A^coii at mid and 
forward/backward rapidity can be projected to Cu+Cu 
and Au+Au using the corresponding parton distribution 
functions for Cu and Au [24|. The results for method 1 
are shown in Fig. |4]as cold nuclear matter baseline Raa 
curves calculated from the best fit values of / (solid lines) 
and the one standard deviation uncertainty in / (dashed 
lines) . The method 2 heavy ion calculations are similar to 
those from method 1, leading to very similar conclusions, 
and are not shown in Fig. [4l In Fig. [4fc) the measured 
Raa values for Cu+Cu are shown divided by the method 
1 calculations for Cu+Cu. The Cu+Cu Raa is seen to be 
consistent with the cold nuclear matter projection within 
about 15% uncertainties up to A^part ~ 50. Given the 
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uncertainty in the cold nuclear matter reference at larger 
-^part values, we can not currently draw any strong con- 
clusions there. However PHENIX completed in February 
2008 a second d+Au run, with approximately 30 times 
the statistics of the first d+Au run in 2003. With the new 
reference d+Au data, we expect to be able to identify if 
and where the measured Cu+Cu -Raa departs from the 
cold nuclear matter baseline. 

In summary, we present high statistics J/tp data from 
Cu+Cu collisions at RHIC, providing for the first time 
detailed information on Raa and (p^) for TVpart < 100. 
The RMS values of the rapidity distributions at all cen- 
tralities are consistent with that for p + p, and the mea- 
sured for pt < 5 GeV/c is nearly independent of 
centrality and rapidity. At similar values of iVpart, ^aa 
and (pt) found to agree within errors for Cu-f Cu 
and Au-|-Au collisions. Cold nuclear matter calculations 
based on ad hoc fits to d+Au data reproduce the periph- 
eral Cu+Cu data well up to iVpart ~ 50, corresponding to 
^Bjorken ^ - 1.5 GeV/fmVc 0, where f-Bjorken is the 
Bjorken energy density and r is the formation time. For 
an estimate of the thermalized energy density, hydrody- 
namical models give thermalization times in the range of 
0.6 fm/c to 1.0 fm/c 0, which implies that cold nuclear 
matter effects dominate J/ip production up to thermal- 
ized energy densities of 1.5 to 2.5 GeV/fm'^. 
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